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Abstract. Observational and theoretical studies point to microquasars (MQs) as possible counterparts of a 
significant fraction of the unidentified gamma-ray sources detected so far. At present, a proper scenario to explain 
the emission beyond soft X-rays from these objects is not known, nor what the precise connection is between 
the radio and the high-energy radiation. We develop a new model where the MQ jet is dynamically dominated 
by cold protons and radiatively dominated by relativistic leptons. The matter content and power of the jet are 
both related with the accretion process. The magnetic field is assumed to be close to equipartition, although it 
is attached to and dominated by the jet matter. For the relativistic particles in the jet, their maximum energy 
depends on both the acceleration efficiency and the energy losses. The model takes into account the interaction of 
the relativistic jet particles with the magnetic field and all the photon and matter fields. Such interaction produces 
significant amounts of radiation from radio to very high energies through synchrotron, relativistic Bremsstrahlung, 
and inverse Compton (IC) processes. Variability of the emission produced by changes in the accretion process 
(e.g. via orbital eccentricity) is also expected. The effects of the gamma-ray absorption by the external photon 
fields on the gamma-ray spectrum have been taken into account, revealing clear spectral features that might be 
observed. This model is consistent to the accretion scenario, energy conservation laws, and current observational 
knowledge, and can provide deeper physical information of the source when tested against multiwavelength data. 

Key words. X-rays: binaries - stars: winds, outflows - gamma-rays: observations - gamma-rays: theory 



1. Introduction The detection of extended non-thermal radio emission 

\ provided clear evidence for the presence of relativistic 

CO , Microquasars are X-ray binary systems (XRBs) with rel- leptong ^ ^ q{ MQgj although it was not consid _ 

atdvistic bipolar outflows or jets (Mirabel & Rodriguez ered ^ gener&1 ^ ^ significantly at x _ 
1999 ). These extended structures have been observed in rayg Qr bey(md Paredes ^ ^ ^ pr0 posed the micro- 
galactic objects atradio wavelengths from theseventies quagar Lg 503Q ag the counterpart of the EGRET source 
(SS 433; Spencer|M3 Hjellming & Johnston [T98TJ). The 3EQ J182 4_ 1514 (Har tman et al. 1999). In their scenario, 
inner region of the disk emits thermally at soft X-rays, the jet relativistic electrons scatter the photons of the mas- 

losing accretion kinetic energy through viscosity-related , ni • , at. a • j. 

° aj ° . sive stellar companion, snowing that microquasar jets are 

phenomena. Additionally, there seems to be evidence sup- .-, , r t-, , , , , . , . , , 

1 1 possible sources ol gamma-rays, further statistical and 

porting the existence of a hot relativistic plasma around theoretical studies showed that microquasars could be be- 

the compact object (called 'corona' hereafter), which is hind gome of the unidentified gamma _ ray sources in the 

radiatively significant during the so-called low-hard state Gakxy (Kaufman Bernad6 et ^gjjjjg Romero et al.TOl 

of the source. Corona electrons would Comptonize soft Bosch _ Ramon et al l flMSalbOOSbl Observational evidence 

disk photons to higher energies (McClintock & Remillard c ■ , , • , ., , c , , , • r 

———T t> t> v j- jets as high energy emitters came from the detection ol 

2004). Transport of angular momentum and kinetic energy X ray ext e nd ed emission (e.g., in SS 433, Marshall et al. 
linked to a specific inner disk magnetic field configuration ^ MigUari ^ &1 2002] XTE J1550 _ 564j Corbel et a l. 
could lead to the forma tion of a jet (MeierEE^. m and H 1743.322, Corbel et al-EHE). The recent de- 
Send offprint requests to: V. Bosch-Ramon tection b ^ the ground-based Cherenkov telescope HESS 
e-mail: vbosch@am.ub.es °f TeV emission coming from the microquasar LS 5039 
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(Aharonian et al. 2005) largely confirms the association 
proposed by Paredes et al. ©000), and it is strong evi- 
dence that microquasars are sources of very high-energy 
gamma-rays, their jets being the best candidates to gen- 
erate such emission. 

For the modeling of gamma-ray emission from jets 
of MQs, there have been two types of approach. One 
considers that hadrons lead radiative processes at GeV- 
TeV gamma-rays and beyond (Romero et al. 1217)51 
Romero & Orellana|[2M5J, producing detectable amounts 
of neutrinos (Torres et al. I2004[) . and leaving electrons 
as possible significant emitters only at lower energies. 
These are the so-called hadronic models. The other ap- 
proach extends the energy of leptons from synchrotron 
radio emitting energies exploring inverse Compton and/or 
synchrotron emission in the jets (i.e. Atoyan & Aharonian 
nM31 Markoff et al. I2UUT1 Georganopoulos et al. I3UU21 
Bosch- Ramon et al. 2005a|. These are the so-called lep- 
tonic models. All these models developed so far are impor- 
tant to investigate to what extent each mechanism of emis- 
sion would be relevant under different circumstances, and 
what level of physical detail is required for realistic model- 
ing with the available observational data. Nevertheless, a 
comprehensive MQ jet model attempting to explain emis- 
sion properties in the whole range of spectral frequencies, 
in accordance with the energy and matter constraints im- 
posed by accretion and the conservation energy law at the 
microscopic level, is still lacking. 

In this paper, we investigate persistent MQ compact 
jets to give mult i wavelength and variability predictions 
consistent with the MQ scenario as a whole. We use the 
term 'compact jet' or 'jet' referring to the type of outflows 
thought to be present during the low-hard state (Fender 
et al. l2003aTl . The extended radio lobes, which are also ob- 
served in some MQs (e.g. IE 1740.7-2942; Mirabel et al. 
1992) , and the blobs ejected during state transitions (e.g. 
GRS 1915+105, Mirabel et al. IT39"S| are not considered 
as compact jets. 'Consistent' means here to develop the 
model taking into account the total amount of matter 
available for accretion, the pressure relationship between 
compact jets and their environments, the standard mod- 
els for accretion and jet ejection, the mechanism for par- 
ticle acceleration, the pair creation and annihilation rates 
and the law of microscopic energy conservation. Semi- 
analytical calculations have been implemented to com- 
pute all the significant emission and absorption mecha- 
nisms that take place in the jet: synchrotron, relativistic 
Brcmsstrahlung with internal and external matter fields, 
inverse Compton with internal and external photon fields, 
and creation and annihilation of pairs. In a forthcoming 
paper, we will present an application of the present model 
to the microquasar LS 5039 (Paredes et al.f2005). 

In Section [3J the general picture of the model is pre- 
sented; in Section [31 the details of the jet model are 
given; in Section the considered radiative processes 
are explained. The resulting spectral energy distributions 
(SEDs) for several relevant situations are shown and dis- 
cussed in Section 03 as are the variability properties of the 



model. Other questions implied by our model are discussed 
in Section [51 and all the treated issues arc summarized in 
Section □ 

2. General picture 

The MQ scenario considered here consists of a binary sys- 
tem formed by a star, either of low or high mass, and a 
compact object, either a black hole or a neutron star. At 
this stage, the nature of the compact object is not rele- 
vant. The stellar companion feeds the accretion process 
onto the compact object. Part of the energy associated 
with the accreted matter is radiated in an accretion disk, 
and part is converted to kinetic and magnetic energy of the 
accretion flow under the effects of the compact object po- 
tential well. In the low-hard state, when the accretion rate 
is moderately low, the accretion disk is well described by a 
standard optically thick and geometrically thin Keplerian 
disk (Shakura & Sunyaev 1973) up to a certain transition 
radius, assumed here to be i?disk ~ 50i?s c h- At distances 
to the compact object smaller than the transition radius, 
we assume the existence of a hot plasma (the corona), 
whose relativistic electrons scatter soft disk photons from 
the geometrically thin disk. The properties of this inner 
region are considered to be suitable for jet ejection phe- 
nomena to take place (Meier I2003|l . The energetics of the 
jet is assumed to be dominated by accretion, and further 
energy sources like compact object rotation have been ne- 
glected (an approximation that seems to be reasonable, 
see Hujeirat 2004 ). A more extended discussion on jet for- 
mation and energetics is present in Sect. l3~T1 

If the jet is formed by accreted matter, it will con- 
tain protons and electrons, as well as a magnetic field (B) 
associated with the plasma. Our assumption is that the 
matter kinetic luminosity is higher than the magnetic lu- 
minosity (or total magnetic energy crossing a jet cross- 
section per time unit) in the jet regions we are concerned 
with, although the magnetic field can be still significant 
once the jet is formed, since the ejection mechanism is 
likely to be magneto-hydrodynamical. For simplicity, the 
jet is supposed to be perpendicular to the orbital plane. 
An important fact is that variability can be easily repro- 
duced when the accretion rate is not constant due to, e.g., 
orbital eccentricity. If protons are relativistic and a con- 
fining mechanism is absent, a mildly relativistic jet (see 
Sect. 13. 2J) will expand relativistically, appearing not colli- 
mated. Therefore, cold protons instead of relativistic pro- 
tons should dominate the jet pressure, with free expan- 
sion speeds smaller than in a full relativistic jet, allow- 
ing for a small jet opening angle. The shock acceleration 
mechanism, assumed at this stage to take place all along 
the jet due to velocity variations in the ejected matter 
(ReesEZHl, can effectively accelerate particles, but must 
such particles diffuse through the shock. This condition 
is fulfilled only by a small fraction of protons that, as- 
suming they follow a thermal distribution, populate the 
highest energy tail. The condition that particles must dif- 
fuse through the shock to be effectively accelerated im- 
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poses a more restrictive condition on the electrons since, 
in general, they have a diffusion mean free path much 
smaller than protons (i.e. lower temperatures), it being 
unlikely that shock acceleration will operate significantly 
on any fraction of the electron thermal distribution that 
also forms the jet. Therefore, we assume that an unspeci- 
fied injection mechanism, like some kind of plasma insta- 
bility, operates on the jet electrons heating them enough 
so they diffuse through the shock (see Sect. 13.2(1 . 

In this scenario, radiative processes that take place 
in the jet can produce significant emission in the whole 
spectrum (see Sects. 0] and EJ. Although there could still 
be enough protons to be significant from the radiative 
point of view, treatments on them can be found elsewhere 
(Romero et al. 120031 |T005 ) , and we will focus on the lep- 
tonic component only 1 . Synchrotron emission can be im- 
portant from radio to soft gamma-rays, if the magnetic 
field is close to equipartition with matter. Moreover, in the 
regions close to the jet ejection point, magnetic and rela- 
tivistic particle densities are high enough for synchrotron 
self-Compton process (SSC) to be considered. Moreover, 
other self-Compton processes can occur, although SSC will 
be dominant in most cases. IC interactions between jet 
electrons and external photon fields (external Compton or 
EC) can also be significant. If the companion star is mas- 
sive, stellar IC scattering can be dominant at the highest 
energies. EC scattering with disk and corona photon fields 
could be non-negligible as well, and relevant for a low- 
mass system (see, nonetheless, Grenier et al. I2005fl . All 
IC fields appear to dominate emission from soft gamma- 
rays to very high-energy gamma-rays (VHE), although a 
relativistic Bremsstrahlung component can also reach very 
high energies, generally being a minor component. Besides 
radiative processes, absorption processes are to be taken 
into account regarding gamma-ray radiation created in the 
jet. Absorption of photons by pair creation can be signif- 
icant at the base of the jet, and also significant within 
the binary system if a massive and very luminous star is 
present. Moreover, secondaries created by electromagnetic 
cascades should be considered, since their contribution is 
not negligible. In this work, we provide a rough estimate 
of these effects on the observed spectrum. It is discussed 
in Sect. 01 and Sect. El For clarity, in Fig. fj] we show a 
rough sketch of the microquasar scenario. 

Four different cases are explored in this paper. First, 
a high-mass and a low-mass XRB presenting a jet with 
a high particle acceleration efficiency are studied in case 
A and B, respectively. For these two particular scenar- 
ios, the accretion disk and the corona have been taken 
to be faint. Here, the effect of the star on the predicted 
spectra is investigated. In case C, the model is applied to 



1 We note that the total luminosity radiated by relativistic 
protons through proton-proton interaction should be close to 
that calculated for relativistic Bremsstrahlung due to similar 
cross-sections and target densities, assuming the same energy 
distribution for both relativistic protons and relativistic elec- 
trons 



observer 




Fig. 1. Picture of the microquasar scenario. The main fea- 
tures considered by the model are listed in their associ- 
ated region: accretion and jet matter rate, magnetic field, 
particle acceleration, energy losses, leptonic emission, pair 
creation and secondary emission. The plot is not to scale. 



a high-mass system considering luminous accretion disk 
and corona with a low particle acceleration efficiency. For 
these three cases the velocity is fixed to a mildly rela- 
tivistic value. This provides restrictions to the fraction of 
accreted matter that is ejected, since the energy budget is 
limited and some part is radiated during the accretion pro- 
cess. Also, the different acceleration efficiency has a strong 
effect on the spectrum at the higher energies. In case D, a 
jet pointing towards the observer without a speed restric- 
tion and with a faint disk and corona is investigated. A 
high particle acceleration efficiency has also been consid- 
ered for this particular scenario, which shows that the jet 
could be a very strong emitter of X-rays and gamma-rays 
under suitable conditions, appearing almost as an ultra lu- 
minous X-ray source (ULXs) if located at large distances. 
We recall that it is supposed that all these systems are in 
a low-hard-like state, when compact jets are thought to be 
present with rather stable characteristics (e.g. Gallo et al. 

Ennii- 

3. The model: formation and properties of the jet 

We describe in this section the simple scenario we adopt 
concerning the formation, collimation and other properties 
of the jet but the radiative ones, which are given in Sect.0J 
The fixed parameters of the model, their representative 
symbols, description and values are provided in Tabled 
being discussed in the text. For the stellar mass loss rate 
(to w ), the jet size, and the orbital parameters, we have 
used typical MQ values. The parameter values that vary 
for the four specific scenarios A, B, C and D are given 
below the fixed ones in Tabled 
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Table 1. Model parameters 



Parameter: description [units] 


values 








e: eccentricity 


0.3 








a: orbital semi-major axis [-R©] 


45 








m w : stellar mass loss rate [Mq yr _1 ] 


10" 6 








kTdisk'- disk inner part temperature [keV] 


0.1 








Pcoi'- corona photon index 


1.6 








^idiak: disk inner radius [i?sch] 


50 








ry. launching radius [J?sch] 


4 








zq'. jet initial point in the compact object RF [-/^sch] 


ou 








X'- jet semi-opening angle tangent 


0.1 








q: equipartition parameter 


0.1 








max. ratio hot to cold lepton number 


0.001 








q acc : fraction of the Edington accretion rate 


0.05 








p: electron power-law index 


2.2 










A 


B 


C 


D 


M x : compact object mass [Mq] 


3 


3 


3 


15 


R+: stellar radius [Rq] 


15 


1 


10 


10 


M*: stellar mass [Mq] 


30 


1 


20 


20 


L*: stellar bolometric luminosity [erg s~ ] 


10 39 


10 33 


10 38 


10 38 


T*: stellar surface temperature [K] 


4 x 10 4 


6 x 10 3 


3 x 10 4 


3 x 10 4 


ft: jet-accretion rate parameter 


0.05 


0.05 


0.02 


0.01 


f: shock energy dissipation efficiency 


0.25 


0.25 


0.05 


0.5 


6: jet viewing angle [°] 


45 


45 


45 


1 


r\: acceleration efficiency 


0.1 


0.1 


0.0001 


0.1 


Odisk: disk radiative efficiency 


0.025 


0.025 


0.25 


0.025 


a COI : corona radiative efficiency 


0.005 


0.005 


0.05 


0.005 



3.1. Jet formation 

3.1.1. Stellar matter accretion 

The accretion rate adopted in this work has been esti- 
mated assuming that the system is accreting at 5% of 
the Eddington value, thought to be typical for accret- 
ing XRBs in the low- hard state (Esin et al. I1997|l . For 
the accretion luminosity to accretion rate ratio, we have 
adopted 0.05c 2 , which for a compact object of 3 Mq yields 
an accretion rate of 6 x 10 -9 M Q yr -1 . Assuming a spher- 
ical Bondi-Hoyle accretion model (Bondi 1952, see also 
Reig et al. 2003), in the case of high-mass MQs, with 
typical stellar mass loss rates of about 10~ 6 M Q yr _1 
and high wind velocities of few times 10 s cm s _1 inferred 
from spectroscopic observations, the mentioned moderate 
accretion rate is hard to achieve. It likely implies that 
some anisotropy in the wind properties occurs in the di- 
rection towards the compact object. We have assumed a 
low wind velocity in the direction of the compact object 
that can be estimated from the stellar mass loss and ac- 
cretion rates quoted above (see also Table obtaining 
values of Sa 10 8 cm s _1 . Moreover, the adoption of an ac- 
cretion model allows us to provide a rough estimate of the 
effects of eccentric orbital motion on emission variability 
as presented in Sect. [3] normalizing the accretion rate to 
the quoted value at phase 0. We adopt the convention here 
that phase is the periastron passage and when the com- 



pact object is in opposition to the observer. We note that 
periastron passage does not correspond to the accretion 
peak, due to effects concerning the composition of veloci- 
ties for both the compact object and the wind. In the case 
of low-mass MQs, orbital variability is not explored here 
because these systems present circular orbits. 

3.1.2. Ejection velocity of the jet 

The jet velocity is estimated taking into account the 
amount of available kinetic luminosity that can be ex- 
tracted from the accretion at the launching radius n (see 
Tableland Sect.[|)|). n is not the distance at which the jet 
is formed, but a characteristic radius where ejected mat- 
ter gets extra kinetic energy from the accretion reservoir. 
The knowledge of such a quantity could help us to under- 
stand where ejection originates, although ejection is likely 
to be a spatially extended phenomenon. It is necessary to 
establish the matter rates of both the advected and the 
ejected matter components that, in relation to the total 
accretion rate, follow the formula: 

m acc = 2m jot + m a dv (1) 

The factor 2 is due to the existence of a jet and its coun- 
terjet. We introduce the parameter k through the rela- 
tionship: TOjct = Km acc . To assign a certain amount of 
extra kinetic luminosity to the ejected matter, we have 



Bosch- Ramon, Romero & Paredes: Broad-band emission from microquasars 



5 



taken into account the energy dissipated in the disk and 
the corona in the form of radiation, as well as the energy 
borne by the advected matter after transferring part of its 
kinetic energy to the jet. This remaining advected kinetic 
luminosity has been assumed to be associated with the 
Keplerian velocity at the launching radius, which is: 



T , x -RSch . 2 
^kadvinj ~ — : "T-advC . 

4ri 



(2) 



The previous considerations give a first order estimate of 
the injection jet velocity. The equation we obtain is: 

jet d~ L]/_ ac j v + idisk + L COT . (3) 

jet is the jet kinetic luminosity at the launching radius. 
This accounts for the kinetic luminosity required to carry 
the jet matter outside the potential well (e.g. it could be 
in the form of magnetic luminosity in those regions, be- 
fore the jet is completely formed at the assumed distance 
of 50 i?sch) plus the kinetic luminosity of the jet after 
ejection: 

GM t 2m j6t 



^kjet(n) 



n 



+ (r jct -l)(2m jet )c 2 



(4) 



(Tj e t — l)(2rhj ct )c 2 , the jet matter kinetic luminosity, ac- 
counts also for the magnetic field and relativistic particle 
luminosity since, once formed, the jet is assumed to be cold 
matter dominated. Ldisk is the disk radiated luminosity, 
taken to be a few % of L acc , and L COI is the corona radi- 
ated luminosity, about 1% or less of L acc (the efficiencies 
that we have adopted are similar to those found in the lit- 
erature; McClintock & RemillardG004): Ldisk = adisk-^acc 
and L cor = a C or^acc- The final expression for the Tj Ct is: 

1 I -Rfieh . -^disk + L c , 



r jet (n) = i + 



2m 



jet 



-RSch . 
— ; — "T-adv 
4ri 



(5) 



We note that the energy dissipated in the shocks 
formed in the jet is extracted from the corresponding 
jet kinetic luminosity after ejection, and the velocity at 
infinity has an associated kinetic luminosity such that 
Lk oo < -^kesc- Therefore, the ejection jet Lorentz factor 
will be larger than its observable value, r^, although both 
will be related by the amount of energy dissipated in the 
jet (see below). This approach is classical, just to provide 
a zeroth order estimate. 

In Sect. [5] the computed SEDs for four different cases 
are presented. Except in the case of D, Tj ct is fixed to 
1.5. This will allow us to compare, through Eq. JSJl, the 
amount of matter carried by the jet between the case 
when disk/corona emission is weak and the case when it 
is strong. This will imply that, for the same ejection ve- 
locity, the jet in the former case can be heavier than in 
the latter one (see the corresponding n values in Table QJ. 

3.1.3. Magnetic field 

Although we do not consider any particular theory of jet 
ejection, it is supposed that the mechanism is a magneto- 
hydrodynamic one. If the jet is ejected by converting mag- 
netic energy to matter kinetic energy, it seems likely that 



B close to the compact object must be beyond equiparti- 
tion with jet matter. We assume that B goes down as it 
transfers energy to the jet matter, accelerating it. Finally, 
B becomes dynamically dominated by jet matter (in the 
jet reference frame (RF)). It is assumed that, when the 
jet is already formed, the magnetic field is entangled with 
matter and approximately turbulent. Both matter and B 
follow adiabatic evolution when moving along the jet, with 
the energy density oc 1/z 2 for a conical jet, where z is the 
distance to the compact object. 

The transition from a magnetic field dominated jet to a 
matter dominated jet should not be discontinuous; B must 
reach values below equipartition in relation to matter, but 
not by too much, since their energy densities should evolve 
in a similar way once the former is attached to the latter. 
A more extended discussion of this issue concerning ex- 
tragalactic jets can be found in Sikora et al. Q2005[l . In 
our work, the magnetic field in the jet reference frame 
at different distances from the compact object has been 
calculated as follows: 



B(z) = y/g 87re p , 



(6) 



where, for a cold proton dominated jet, the jet matter 
energy density is: 



e P 0) 



TOjet 



^jet^etmp 



<E pk >- 



2ttz 



(7) 



We introduce R; 



jet 



XZ, and < E p k >, which is the 
mean cold proton kinetic energy, taken to be the classi- 
cal kinetic energy of protons with velocity equal to the 
expansion velocity (V exp = x^jct)- This would correspond 
approximately to the sound speed of the plasma in the jet 
RF. 



3.2. Jet properties 
3.2.1. Confinement 

Compact jets in microquasars appear to be in general 
mildly relativistic (for LS 5039, see Paredes et al. 2002, for 
LS I +61 303, see Massi et al. 120041 in general, see Gallo 
et al. I2003|l . This implies that these jets cannot be colli- 
mated by relativistic motion, and an external or internal 
factor must collimate them. External medium collimation 
operates when its pressure is similar to or larger than the 
jet lateral pressure, both taken as a non relativistic ideal 
gas (the jet expansion velocity considered here, as well 
as the environment gas, is not relativistic). Environment 
pressure can be exerted by expelled stellar matter (e.g. 
stellar wind) and the interstellar medium (ISM). In Fig. [2 
we show the jet pressure along the jet, the stellar wind 
pressure and the ISM pressure. The jet specific parameter 
values relevant here are those corresponding to the case 
A. For the ISM, we have adopted a density of 10 cm~ 3 
and a temperature of 10 3 K, assuming that it is heated 
by the massive star. For this particular case, the wind ve- 
locity has been fixed to 10 s cm s _1 . As is seen in Fig. |3 
the pressure of the wind is not significantly above the jet 
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pressure anywhere along the jet. There is a clear change 
of behavior of the wind pressure on the jet when it turns 
from kinetic to thermal pressure, at distances similar to 
the binary system size. We note that ISM effects could be 
non-negligible at very large scales. We have not studied 
the interaction of the jet with the ISM, although a treat- 
ment of this kind can be found elsewhere, e.g., in Heinz 
( 2002 ). Further sources of external influence could be disk 
radiation and/or disk winds. However, the former does not 
appear to be effective at middle and large jet scales (see, 
e.g. Fukue et al. [2001 ) , and the latter, if effective beyond 
disk scale distances, would need to be as compact as the 
jet itself, and the problem of collimation is transferred to 
this compact disk wind. 

Internal collimation could be achieved by a special 
magnetic field configuration. Nevertheless, since the mag- 
netic field energy density is considered to be lower than 
the matter energy density (see Sect. 13.1} . this mechanism 
does not seem to be relevant here. In our context, the jet 
expands freely at roughly its sound speed, and this type 
of jet can only be collimated if the plasma sound speed is 
low, implying a cold matter dominated jet. For instance, 
the jet of SS 433 is particularly heavy among MQ com- 
pact jets and seems to be cold matter dominated, since 
there is significant thermal X-ray radiation coming from 
the jet (see, i.e., Fender et al. I2003b| . Jet opening angles 
are typically of a few degrees (for LS 5039, see Paredes 
et al. l2"UU2l for SS 433, see Marshall et al. This 
means that the parameter \ must be similar to or less than 
0.1 (see Table nj. For mildly relativistic jets, if their cold 
components are in thermal equilibrium, the temperature 
associated with V cxp will be similar to that of the inner 
regions of the accretion disk or corona (several 100 keV). 

In order to obtain some upper limits for the amount 
of accelerated particles, we have computed both the cold 
and the hot proton pressure along the jet. For hot par- 
ticles, we have estimated their pressure as corresponding 
to a relativistic ideal gas of protons following a power-law 
distribution. It is found that the pressure of the hot com- 
ponent is similar to that of the cold component for a ratio 
of hot particle to cold particle number of about 1 /1000. We 
have adopted this value for £, the parameter that gives the 
maximum possible ratio of hot to cold lepton number (we 
assume the existence of one cold lepton per cold proton). 
We note that the energy density of relativistic leptons is 
below the magnetic one all along the jet. This fact could 
imply a better confinement for such particles, considering 
that the gyro-radius becomes smaller as the magnetic field 
increases. 

3.2.2. Relativistic particles 

Internal shocks due to different velocities of the jet plasma 
(ReeslMHl Spada et aL ESUUl Yuan et al. l2Tjtj5l> can dissi- 
pate bulk kinetic energy, converting it into random kinetic 
energy of accelerated particles. The injected particle spec- 
trum is kept all along the jet up to a certain maximum 
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Fig. 2. Pressure values along the jet for the jet matter 
(solid line), stellar wind (dashed line) and the ISM (dotted 
line). To compute this, k = 0.05 and a ISM with a density 
of 10 cm -3 and a temperature of 10 3 K have been adopted. 



energy, which varies in accordance with the balance of 
energy gains and losses. In addition, the assumption that 
the jet is cold matter dominated stringently constrains £ to 
the quoted value. Effects of cooling on the power-law spec- 
trum cannot be considered at this stage due to the uncer- 
tainties concerning the injection/acceleration processes, 
and its shape index is therefore fixed. Accelerated leptons 
emitting synchrotron radiation can produce jet emission 
similar to that observed at radio wavelengths in LS 5039 
(Paredes et al. 2000), the microquasar detected by HESS. 

Under the conditions assumed here, the first order 
Fermi mechanism cannot accelerate thermal electrons 
fBcll ll97%j) efficiently. Therefore, a still unspecified mech- 
anism that can accelerate thermal electrons must operate 
up to the minimum Fermi acceleration energy, 7 m i n . To 
determine 7 m i n , we have assumed that the accelerated lep- 
tons have a mean free path within the shock similar to the 
mean free path of the particles that form the shock, i.e. 
the jet cold protons. The injection mechanism remains un- 
specified (e.g. magnetic reconnection that heats electrons, 
some type of proton-electron temperature coupling, injec- 
tion of relativistic pairs, etc), and the injection rate will 
depend on acceleration efficiency constraints. The Fermi 
mechanism provides the energy of radiating leptons, and 
not this unknown injection mechanism. Thus, the radia- 
tive properties do not rely on a mysterious process of par- 
ticle injection but on the well-known first order Fermi the- 
ory of particle acceleration (a similar approach to hadron 
acceleration is performed in Mastichiadis & Kirk|l995). 

The relativistic lepton energy distribution in the RF 
of the jet is assumed to follow a power-law: 

"c(7,*)= 7~ P , (8) 
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where iVj n j is the injection normalization parameter for the 
energy distribution of relativistic leptons, which is taken 
to vary as a function of the acceleration mechanism effi- 
ciency and jet conditions (see below). The injected lepton 
rate Qmj (i- e - the total number of leptons crossing a jet 
section at z per second) is associated with N^, and can 
be obtained from the latter integrating N- m jj~ p from 7 m ; n 
to 7 m ax- P is taken to be 2.2, a reasonable value that can 
be derived from optically thin radio spectra observed in 
some MQs. Since the matter density decays as 1/z 2 in a 
conical jet, we observe the z-dependence in Eq. (JSJ. The 
acceleration process is assumed to keep the same energy 
distribution for the relativistic leptons along the jet, al- 
though the maximum energy of the accelerated particles 
depends on energy loss conditions and the size of the ac- 
celerator, taken here to be the jet width. We will neglect 
at this stage the effects of escaping particles on the par- 
ticle energy distribution (for a detailed treatment on this 
see Atoyan & Aharonian 1999). The power-law is simply 
cut at a certain 7 max , which is computed as explained 
below. Nevertheless, since the accelerator/jet size limits 
the acceleration efficiency, escape losses are taken into ac- 
count as cooling terms. We note that a low energy particle 
population could be present in the jet just below the ac- 
celeration injection energy, that can have Lorentz factors 
^100. These particles cannot radiate significantly, due to 
their long radiative timescales and because of adiabatic 
losses, at least while they are confined within the jet. 

3.2.3. Particle acceleration and jet deceleration 

We have adopted Fermi first order acceleration theory to 
calculate the maximum energy of accelerated particles in 
the jet (see, e.g., Biermann & Sttritmatter ll987l Protheroe 
1999J. We set the acceleration rate equal to the rate at 
which particles lose their energy. Expressed in terms of 
the Lorentz factor: 



Tgain — 7loss ■ 

The energy gain rate can be calculated from: 
rjq c Bc 



7gain 



TO C C Z 



(9) 



(10) 



where r\ can have different values depending on the shock 
conditions 2 . For simplicity and since the specific shock 
conditions are not known, r\ is treated here as a free pa- 
rameter in the range 10 _4 -10 _1 . r\ is likely a function of 
z, although at this stage is set to be constant. We recall 
that B goes down by 1/z. The energy loss rate can be es- 
timated adding the contribution of the different types of 
energy losses. On the one hand, there are adiabatic losses: 



7adiab 



ill 



2 It depends on the angle between the magnetic field lines 
and the perpendicular direction to the shock surface, and on 
the shock speed in the plasma frame, as well as the diffusion 
coefficient in the shock region, rj is typically between ~ 10~ 4 — 
KT 1 (Protheroe ESSU 



since relativistic leptons are exerting work against the jet 
confining walls, whose expansion is led by cold protons. 
On the other hand, there are radiative losses due to syn- 
chrotron, IC and relativistic Bremsstrahlung processes. 
Hence, we obtain: 



7io 



Trad + 7adiab- 



(12) 



The expressions for the particle energy loss rates for each 
radiative mechanism can be found, e.g., in Blumenthal 
& Gould 1)1970(1 . From these expressions together with 
Eq. (O, one can finally obtain 7 max and its evolution with 
z. 

There is a limit for the total amount of energy avail- 
able for shock acceleration, which is likely related to the 
shock efficiency to dissipate energy via heating jet leptons. 
The upper limit for shock energy dissipation efficiency is 
that 100% of the kinetic energy is dissipated. It seems un- 
likely that a so efficient process takes place but, still, an 
important amount of the energy is released in the form 
of radiation, as it appears to happen for extragalactic jets 
and some galactic jets. We have adopted the criterion that, 
at most, the shock could dissipate up to some fraction of 
the flow average kinetic energy accelerating particles in 
the whole jet, or £ = 0.05, 0.25, 0.5 (see Table QJ. It is not 
stated that a 5, 25 or 50% of the whole jet kinetic energy 
goes to the relativistic particles. The shock could dissi- 
pate enough energy to accelerate the maximum number 
of relativistic leptons in the jet (given by (), still being 
below the efficiency limit. This implies that for most of 
the jet the energy dissipation rate can be below this con- 
straint. If a large fraction of the jet kinetic energy goes to 
heat particles, the shock dynamics are affected. We do not 
include such effects in our calculations, treating the accel- 
eration process as the test particle case. We do not include 
the effects of the shocks in the local conditions either, but 
one can consider them as averaged ones (for more precise 
calculations, properties of the shocks like velocity or com- 
pression ratio would be required, which is beyond present 
knowledge and it is absorbed by other parameters like rj 
and g). £ might be interpreted as referring to a feedback 
effect on shock dynamics produced by accelerated parti- 
cles. 

The next constraint has been imposed: 



dL Ais (z) dL\ oss {z) 



(13) 



dz dz 

where dL\ QSS (z)/dz is the energy lost by length unit at 
different z through radiative and adiabatic losses and has 
the form: 



dLi oss (z) 
dz 



irRf ct n c ('y,z)ji oss dj, 



(14) 



provided that all the quantities are already in the com- 
pact object RF. dLdi s (z)/dz is the maximum dissipated 
luminosity per length unit, which is taken to evolve with 
z like 7 gain : 



dLdi S (z) Cdi f 
dz z 



(15) 
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The normalization constant Cdis can be obtained by in- 
tegrating the total amount of shock energy available per 
length unit along the jet, from Zq to z max , and equating it 
to a suitable fraction (£) of the jet kinetic luminosity: 



38.0 



case A & B 
case C 
case D 



M-Zmax/^o) ' 



(16) 



CO 



37.0 



36.5 - 



dLdis(z) / dz has a weak 2 max -dependence. This quantity 
has been taken to be about 0.1 pc, where environmen- 
tal effects on the jet properties could become significant, 
changing the characteristics of the jet, as it has been ar- 
gued in Sect. 13.21 

Equating Eq. |T3J) and Eq. JTSJ (through Eq. (SJ gives 
iVj n j in the compact object RF, although this value can- 
not be higher than that derived from £ either. In Fig. [21 
we show the evolution of Qinj for the four MQ cases stud- 
ied here. Losses are strong enough well within the binary 
system to hold the total number of relativistic leptons in- 
jected per second (i.e. crossing a jet section at z) below the 
maximum allowed rate. As these particles are carried by 
the jet to larger z, losses are lower and the injected particle 
rate reaches its maximum value. For the cases considered 
here, synchrotron losses dominate over SSC, EC and rela- 
tivistic Bremsstrahlung losses, although star IC losses can 
dominate at certain jet regions and/or for slightly smaller 
values of B. Therefore, the Klein-Nishina (KN) effects are 
to be taken into account. 

In Fig. 0J we see the evolution of the maximum energy 
of the relativistic particles with z. Radiative losses limit 
this maximum energy at spatial scales similar to the size 
of the binary system. At middle and large jet scales, the 
acceleration process is limited by the typical size of the 
jet and by the jet magnetic field. Both z-dependences get 
canceled, and the consequence is a constant value for 7 max . 
We show the evolution of 7 max for two particularly inter- 
esting cases (Fig. 0| : A and C. Since the amount of energy 
transferred to relativistic particles can be significant, we 
calculate the effects of this process on the jet Lorentz fac- 
tor. Since different velocities in plasma motion generate 
shocks and the acceleration process, the jet Lorentz factor 
must be understood as an average value. The results for 
cases A, B, C and D are shown in Fig. 03 As is expected, 
the lower the acceleration efficiency (i.e. the transfer of 
energy from the bulk motion to the radiating particles), 
the lower the decrease in the bulk motion Lorentz factor, 
which goes down strongly at the base of the jet and is 
already stabilized at binary system scales. The radiative 
efficiency of the jets for our particular parameter choice 
(for g, £ and 77) is of about 1-10%, not far from estimates 
obtained by different approaches (see, e.g., Fender [2001 
Yuan et al. 2005). With the Lorentz factors given in Fig.[5| 
the initial jet kinetic luminosities are 8.7 x 10 36 erg s -1 , 
8.7 x 10 36 erg s -1 , 3.5 x 10 36 erg s" 1 and 4.9 x 10 37 erg s" 1 
for A, B, C and D respectively. 
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Fig. 3. The evolution of Q in j along the jet for A, B, C 
and D. The free parameter values adopted for these cases 
are shown in Tabled 
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Fig. 4. The same as in Fig. but for j max in the cases A 
and C. 

4. The model: radiative processes in the jet 

4.1. Radiation mechanisms 

We have accounted for synchrotron, relativistic 
Bremsstrahlung and IC emission (in both Thomson 
and KN regimes). Other leptonic radiative processes such 
as free-free transitions, thermal Bremsstrahlung, etc, 
have been neglected since jet emission is dominated by 
the non-thermal processes, although for a completely cold 
and powerful jet these processes would have observable 
effects. We show in Fig.Elan example of the energy lost by 
leptons per volume and time units in the jet at different 
distances from the compact object. As is seen there, the 
dominant type of loss is synchrotron radiation, although 
the IC losses can be significant. Well outside the binary 
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Fig. 5. The same as in Fig. El but for rj ct normalized to 
its initial value in the four cases explored here. 
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Fig. 6. The same as in Fig.[3]but for the energy per volume 
and time units lost by leptons, de\ osa /dt, now only for A. 



system, the dominant losses are the adiabatic ones. 
For certain ranges of the parameter values, stellar IC 
(A-like but with low magnetic field) and corona IC (C) 
could dominate over synchrotron radiation. Relativistic 
Brcmsstrahlung is negligible in general. Since the increase 
of the amount of relativistic particles stops when their 
number reaches a maximum value (recall £ and Fig. |3J) , 
the decrease in the total emissivity becomes more serious 
at distances of about 10 11 cm from the compact object. 

Taking into account the physical conditions along the 
jet, we have computed the SEDs corresponding to the 
different mechanisms mentioned above. To calculate the 
spectrum of the radiation coming from the jet, we have 
used the energy distribution function of the relativistic 
leptons shown in Sect. 13.21 as well as the corresponding 
cross-sections of each process. For synchrotron emission, 
the magnetic field at different z is given in Sect. I3~T1 and 
we adopt the width of the jet as the length that deter- 
mines whether synchrotron emission is either optically 
thin or optically thick. For external Bremsstrahlung, in 
the case of interaction with the stellar wind ions (the 
wind is considered as a completely ionized plasma), we 
take our calculations as an upper limit and assume that 
the target ion density is the wind density (i.e. that all the 
available wind particles diffuse within the jet); for inter- 
nal Bremsstrahlung, the target ion density is the proton 
density of the jet, directly derived from mj e t since the jet 
is cold matter dominated. For EC emission, the target 
photon densities from the star, the disk and the corona 
are those described in Sect. 14.21 and in the case of inter- 
nal Compton emission (basically, SSC), the target photon 
density has been calculated previously (i.e. synchrotron). 
All the z-dependences have been taken into account, di- 
viding the jet in slices, each with homogeneous properties, 
and the overall jet emission has been integrated over all 
the slices. For further details concerning synchrotron, rel- 



ativistic Bremsstrahlung and IC processes, as well as the 
electron energy loss expressions, we refer to the work by 



Blumenthal & Gould 1)1970(1 . The Doppler boosting effects 
in the observed spectra have been implemented as usual 
(e.g. Dermer & Schlickeiser 2002). 

4.2. External photon fields 

The external seed photon sources considered in the model 
are the star, the disk and the corona. The star and the 
disk have been considered to be gray bodies, normalized 
to their total luminosities. The star photon distribution 
peaks at optical-UV energies, and the disk one peaks at 
100 eV. This disk photon energy appears to be typical in 
the low-hard state when the optically thick disk is trun- 
cated far from the compact object and the disk matter 
does not reach temperatures as high as during more in- 
tense accretion states, when the inner disk radius shrinks. 
For the corona emission, we have assumed that it follows 
a power-law plus an exponential cut-off, peaking around 
100 keV. Since IC interactions are studied first in the jet 
RF, we calculate the total energy densities and the spec- 
tral energy densities in this RF. Since expressions for the 
star, disk and corona photon energy density are originally 
in the compact object RF, they are transformed using a 
relationship found elsewhere (e.g. Dermer & Schlickeiser 
2001: 



r?et(l + /3M) 3 



(17) 



where f2 represents the photon direction, /i is cos $ is 
the angle between the photon direction and the jet axis, 
and eo is the seed photon energy for IC interaction. The 
quantities with primes are in the compact object/observer 
RF, and in the jet RF otherwise. For the star photons we 
assume that they reach the jet, at a particular z', with the 
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same direction, which depends on the orbital phase and 
z' . In the case of disk photons, their direction is taken to 
be coming from behind the jet. To treat mono-directional 
seed photon fields, we have adopted the approach used by 
Dermer et al. I)1992|l 3 . Since the interaction angle between 
jet electrons at different z' and stellar photons affects the 
star IC emission and varies with the orbital phase, we re- 
call the adopted criterion that the systems treated here are 
observed in such a way that at phase or periastron pas- 
sage, the compact object is in opposition to the observer. 
To apply the model to particular objects, the specific rela- 
tionship between observer line of sight and compact object 
phase is needed. For corona photons, the jet base in this 
model is assumed to be located approximately in the ex- 
ternal parts of this region. This source of photons is closer 
to the jet than the disk, and the interaction is significant 
only at the inner region of the jet due to the 1 / z 2 -decrease 
in density of corona photons. Thus, we have adopted the 
assumption that the corona field is roughly isotropic. 

4.3. Pair creation and annihilation within the jet 

We have investigated the effects of pair creation and an- 
nihilation phenomena on the jet emission by calculating 
the gamma-ray opacities (r 77e ') for jet photons of differ- 
ent e' and produced at different z' . To calculate opacities 
as well as the number of created and annihilated pairs, we 
have used the pair creation and annihilation rates given 
in Eqs. 3.7 and 4.6 of Coppi & Blandford flMty . The evo- 
lution of the opacity is shown in Fig. It is seen that 
opacities are higher in the vicinity of the compact object 
at energies between 1 GeV-10 TeV. Gamma-ray absorp- 
tion in the stellar UV photon field is significant for photon 
energies of ~ 100 GeV (for a deeper treatment of gamma- 
ray opacities due to the stellar photon field, see Romero et 
al. 13(7)51 Bottcher & Dermer HHH and Dubus . For 
certain parameter values, the opacity could be significant 
as well at 10 MeV and a few GeV within the corona and 
disk fields (see also Romero et al. |2002 ) . 

The gamma-ray opacity by pair creation inside the jet 
is very high at z' ~ z' Q because the target photon density 
is large and very sensitive to the jet width. If the jet width 
were larger, the internal opacity would be weaker, without 
affecting very much the overall production spectrum. Due 
to these uncertainties, the pair creation due to jet internal 
fields should be studied in a more detailed model in future 
work. 

Concerning annihilation rates inside this cold matter 
dominated jet, for any reasonable set of parameter val- 
ues the luminosity that could be emitted in form of an 
annihilation line is too low to be distinguished from the 
continuum emission. Other models, like the one of Punsly 



3 For the mildly relativistic jets treated here, the disk IC ra- 
diation coming from the counterjet is enhanced by the angular 
dependence of the IC interaction. However, for the adopted pa- 
rameter values, this component is at most similar to the normal 
jet component. 
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Fig. 7. Evolution of r 77e / at different z' . We have adopted 
the parameter values of the case A (see Table 0. 



et al. i|2000[) . where a pure pair plasma is assumed, could 
produce detectable annihilation lines. 

Observable predictions from considering pair creation 
phenomena in our model are presented and discussed 
briefly in Sect. OD and |HJ although we remark that the 
creation of pairs inside the jet could lead to the appear- 
ance of bumps due to the accumulation of particles at 
the energies of pair creation. To introduce such an effect 
properly requires a better knowledge of the particle in- 
jection function, which is beyond the scope of this work. 
Therefore, the high-energy gamma-ray band of the com- 
puted SEDs probably gives good enough flux estimates, 
although slopes could be slightly different as a result of 
all these subtle effects. 

For those pairs that are created within the binary sys- 
tem, but outside the jet, the situation is different from 
that of pairs created inside. Starting with a determinate 
number of relativistic particles in the jet, plus the given 
jet conditions, one can consistently derive the SED of the 
produced radiation in the compact object RF. Thus, the 
spectrum is known, and it allows us to know precisely the 
number of absorbed photons and created pairs within the 
star, the disk and corona photon fields (for previous treat- 
ments of this, see Romero et al. I2002fl . From the previous 
result, it is possible to roughly estimate the radiation that 
is generated by those pairs through IC interaction with 
external source photons. Although it is a rough estimate, 
it is found to be in agreement with more detailed models 
of these processes (Khangulyan & Aharonian 2005 ) . 

5. Application of the model 

The different radiation components produced in the jet 
and the predicted SEDs have been computed for the four 
specific scenarios considered here. The effects of pair cre- 
ation phenomena due to the external photon fields inter- 
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acting with the produced gamma-ray photons are taken 
into account, and the secondary radiation produced by the 
created pairs is estimated. The calculations are performed 
at the periastron passage, when the compact object is in 
opposition to the observer and the interaction angle be- 
tween star photons and jet leptons implies more luminos- 
ity for the star IC component (see Dermer et al. I1992fl . 
showing the importance of such an effect. However, such 
an angle depends on the electron energy, which should 
be taken into account in more detailed models of the IC 
interaction (e.g., Khangulyan & Aharonian[2005). 

The broad-band SEDs for cases A and B are presented 
in Figs. |H1 and [5] respectively. The strong effects on the 
computed SED due to the presence of a massive star can 
be appreciated. The star IC component is very significant, 
partially because of the specific interaction angle between 
seed photons and leptons at phase 0, and also because 
the interaction with stellar photons is more significant at 
large z, when the number of relativistic particles is higher 
(see Fig. than for other photon fields. For A, gamma- 
gamma opacity is very high at VHE. We recall that the 
disk and corona emission have been assumed to radiate 
just a few per cent of the accretion power. As accretion 
does not dissipate a significant fraction of the available 
energy via either disk or corona radiation, the jet can carry 
more energy and matter for the same ejection velocity 
(and the assumptions put forward in Sect. ETTl are valid). 
The acceleration efficiency has been assumed to be high. 

In Fig. EH the broad-band SED of case C is shown. 
We have increased the disk and the corona emission, fix- 
ing the jet velocity. This implies a lighter jet than in the 
two previous cases. Also, we have modified the accelera- 
tion efficiency of the jet for this particular case to a smaller 
value than those previously used. The system harbors a 
massive star, although neither as massive nor as bright as 
in A. In the case of D, whose SED is plotted in Fig. ^2 
the matter content of the jet has been reduced to a smaller 
value, fixing again disk and corona emission and particle 
acceleration efficiency as in A and B. This implies that the 
jet motion will be more relativistic than in previous cases. 
This is unusual in what has been said in this work. Our 
purpose in studying this particular scenario is to show that 
an ultra luminous X-ray source can be reproduced through 
strong beaming and small viewing angles. Since the goal 
in this specific scenario was to obtain very high X-ray 
fluxes, the compact object has been considered a 15 M Q 
black hole, allowing for higher accretion rates since the 
Eddington limit is higher for more massive accreting ob- 
jects. We have also increased the shock maximum energy 
dissipation efficiency up to £ = 0.5. 

The SEDs obtained at high energies resemble roughly 
those obtained applying the model presented in Bosch- 
Ramon et al. (|2D05a). This is because the IC emission 
dominates at high energies and, at least in some of the 
cases, the dominant IC component is the same as in that 
model. However, the physical motivation of the present 
work goes much further than before, with predictions con- 
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Fig. 8. Case A computed SED for the entire spectrum 
as it would be observed. Attenuation of the jet photons 
due to absorption in the external photon fields is taken 
into account, as well as the IC emission of the first gen- 
eration of pairs created within them. Isotropic luminosity 
is assumed. The different IC, relativistic Bremsstrahlung, 
synchrotron and other seed photon fields are shown. For 
the several components, the production SED is shown. 
The corona photon field is also taken into account, but 
its effects on pair creation and subsequent emission are 
overcome by the synchrotron emission. 



cerning radio, variability and jet physical properties that 
could not be provided in the previous model. 

5.1. Spectral properties 

At gamma-rays, a different component is the main one in 
each explored scenario. For A, the star Comptonized pho- 
tons are dominant, reaching 10 34 erg s _1 at 100 MeV, with 
a photon index of about 2, and 5 x 10 32 erg s _1 at 100 GeV. 
At VHE, due to the gamma-gamma absorption, the model 
predicts a soft photon index that hardens at higher ener- 
gies. There could be additional spectral features related 
to the secondaries created in the corona, the disk and the 
star photon field. These features would appear as bumps 
at energies of a few GeV for disk pairs and a few tens of 
GeV for star photons (negligigle in A). Even for negligible 
disk and corona sources, a decrease in the predicted flux 
beyond 50 GeV is unavoidable if it is a relatively close 
and massive binary system. In the case of B, the dom- 
inant components are the corona IC and the SSC ones, 
with luminosities of about 10 33 erg s _1 at 100 MeV with 
a photon index of 2 (beyond the synchrotron component, 
see Pig- EJ) 5 and few 10 32 erg s _1 at 100 GeV. In this case, 
however, the low mass star photon field does not signifi- 
cantly affect the VHE spectrum, and only a small bump 
due to the secondaries created within the disk photon field 
is visible in the SED. Therefore, the photon index beyond 
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Fig. 9. The same as in Fig. [HI but for the case B. The 
small bumps present from beyond 100 MeV come from the 
IC radiation emitted by those leptons generated by pair 
creation in the disk and the stellar photon field. These pair 
components are not made explicit in the plot for clarity. 
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Fig. 10. The same as in Fig. [HI but for the case C with 
a jet with low acceleration efficiency. For this particular 
situation it is possible to see the small bump at 10 MeV 
produced through IC scattering by the pairs created in 
the corona photon field. 



10 GeV gets softer because of the KN effect, although it is 
a smoother steepening than in A. Concerning C, the lumi- 
nosity dominated by the corona IC emission at 100 MeV 
is about 10 33 erg s _1 and the photon index is similar to 
the one predicted in the previous cases. Otherwise, only 
disk and corona gamma-gamma absorption effects are sig- 
nificant, visible in the plots at about 1 GeV and 10 MeV 
respectively (see Fig. llQIi . The star IC component turns to 
dominates beyond a few GeV. For D, the Doppler boost- 
ing affects synchrotron radiation, which reaches 100 MeV 
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Fig. 11. The same as in Fig.[H]but for the case D. 



with luminosities higher than 10 erg s . However, the 
SSC component is not dominant at all because it is very 
sensitive to the jet density, and the jet now is the least 
dense among the four studied cases. Therefore, beyond 
1 GeV the source is dominated by the corona and star IC 
components, with a photon index softer than 2 and some 
absorption produced by the star photon field. Small peaks 
of the IC emission of the secondaries are also present pro- 
duced within the star and the disk photon fields. 

In X-rays, the emission is synchrotron dominated for 
most cases (except for C). This is because the magnetic 
field is below but not far from equipartition. The mat- 
ter energy density at the jet base is so high that the 
magnetic field reaches values around 10 5 G, allowing syn- 
chrotron radiation to dominate up to soft gamma-ray en- 
ergies. Similar results have been obtained by Markoff et al. 
1)20011 I2003[) for the case of shock acceleration limited 
mainly by synchrotron losses. In our model particle accel- 
eration is limited by shock energy dissipation efficiency, jet 
size, adiabatic and all the radiative losses (including KN 
regime for IC losses) . Because of the evolution of 7 ma x that 
rises when z is larger (for reasonable parameter values), 
the synchrotron spectrum changes smoothly at energies 
around 100 keV. Nevertheless, the disk and corona are 
not negligible in general, and in the case of C, disk and 
corona overcome the jet radiation up to 1 MeV. D as an 
X-ray source is extremely bright. This result shows that a 
light and fast jet observed from very small viewing angles 
might turn out to be an ULX. Because of the low proba- 
bility of being observed, it is more likely to detect them in 
other galaxies (for previous works on MQs as ULXs, see, 
e.g., Georganopoulos et al. 120021 Kording et al. l2002|l . 

Below stellar emission energies, synchrotron radiation 
dominates again. If stellar emission were reprocessed by 
absorption and shifted to lower energies, it is likely that 
the far infrared would still be dominated by the en- 
shrouded stellar component. At radio frequencies, there 
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is significant radio emission, with isotropic luminosities of 
about 10 29 erg s _1 at 5 GHz (few mJy at 3 kpc). Our 
model predicts core-dominated emission, strongly self- 
absorbed due to the high efficiency of the inner jet to 
radiate through synchrotron process. The SED in the ra- 
dio band corresponds to a spectral index equal to zero (or 
eL e oc e +1 ), as it would be expected from a conical jet. 
C presents the weakest radio emitting jet, weaker than 
in A and B, and D is the strongest radio source due to 
the Doppler boosting. Further aspects on radio emission 
in our scenario are commented in next section. 

5.2. Variability 

Variability through changes in the stellar mass-loss den- 
sity profile is introduced in the model in a consistent way 
when referring to high-mass microquasars. In Fig. El the 
SEDs of A at phases 0.3 (accretion maximum), (peri- 
astron passage) and 0.7 (accretion minimum) are shown. 
For clarity, we have split the overall spectrum in three 
bands: radio, X-rays and gamma-rays. Emission varies due 
to orbital eccentricity for a spherical slow wind accreting 
system with constant wind velocity. Even for the simple 
model considered here, it is possible to see from the figures 
the complex evolution of radiation. At radio and X-ray 
frequencies, the different fluxes associated with different 
phases are correlated with the accretion rate. Since the 
velocity of the compact object relative to its surround- 
ing medium changes along the orbit, the periastron is 
not directly associated with the highest flux. However, at 
gamma-rays, the particular angle of interaction between 
star photons and jet electrons implies that the stellar IC 
is dominant, this phase being associated with a peak in 
gamma-ray emission. More complex variations in the ac- 
cretion, affecting jet ejection and radiation, could be intro- 
duced through the function reproducing the stellar mass- 
loss density profile. Other timescales of variability, linked 
to disk phenomena, could be introduced through the pa- 
rameter that controls the amount of matter that goes to 
the jet, although it is beyond the scope of this work. Also, 
the radio variability pattern depends on the scale of the 
radio emitting region, which is the inner jet for the pa- 
rameter values adopted here. 

6. Discussion 

We conclude that persistent jets in MQs in the low-hard 
state, despite accreting at relatively low rate, and un- 
der reasonable conditions for the jet matter, energy and 
magnetic field, can radiate with significant luminosities 
from radio to gamma-rays. The model provides predictions 
about the shape of the SEDs and points to MQs as VHE 
sources, as has been recently confirmed by Aharonian et al. 
( 2005J) in the case of LS 5039. It also predicts variability 
at different ranges of energy. The importance of the syn- 
chrotron cooling channel for relativistic particles in the jet 
is high. This points to the fact that, for certain objects, the 
jet could overcome at all wavelengths any other emitting 
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Fig. 12. The case A, three predicted SEDs at radio, X-ray 
and gamma-ray energy bands, corresponding to three dif- 
ferent orbital phases: (periastron passage, dashed line), 
0.3 (highest accretion rate, solid line), and 0.7 (smallest 
accretion rate, dotted line). 

region of the MQ, except the star itself for high mass sys- 
tems. We note that for a relatively weak disk and corona 
(10 34 erg s _1 ), and even with g ^0.01 (one tenth of the 
one adopted here), the spectrum would be jet dominated 
and well described by a power-law at X-ray wavelegnths 
(see also Paredes et al. 12005). An observational feature 
that could determine whether X-rays come from a jet is 
the dependence of X-ray fluxes on the accretion rate, if the 
latter can be estimated. In the scenario presented here, X- 
rays vary with m acc , but it is also indirectly connected to 
B and 7 ma x, the latter also being sensitive to the over- 
all jet conditions and size, and both quantities depend on 
the accretion rate as well. Long exposure multiwavelength 
observations will be required if a relationship between ra- 
diation components of different origin is to be found with 
high accuracy, particularly considering the moderate X- 
ray fluxes of permanent jet sources like LS 5039, which do 
not appear to follow the typical behavior of X-ray binaries 
in the low- hard state (Bosch- Ramon et al. 2005c). An ad- 
ditional observational feature that would favor the jet as 
the origin of the X-rays would be the detection of some 
amount of polarization in this radiation. 

At gamma-rays, our jet model, more detailed and in 
accordance to the current knowledge of MQs than our 
previous works, still shows that MQs could be behind 
some EGRET sources (likely those variable and located 
in the galactic plane), and also predicts a different evo- 
lution of emission at different energies due to IC angular 
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dependence interaction, some strange features like radia- 
tion bumps and depressions in the spectra due to gamma- 
gamma absorption with external photon fields, and a non 
trivial relationship between pair creation and particle in- 
jection within the jet itself. GLAST and AGILE, with 
a sensitivity at least several times better than that of 
EGRET, should be able to detect microquasars even when 
they are low mass systems (case B) and/or they have 
weaker jets (case C). The detection of the microquasar 
LS 5039 by HESS shows that the efficiency of the particle 
acceleration processes should be high. Due to the strong- 
photon absorption beyond ~ 100 GeV, it might be that 
the bulk of the TeV emission came from regions where 
stellar photon density is not significant. 

Although we are using accretion rates which are a 
small fraction of the Eddington luminosity, the wind ve- 
locity is required to be low to power accretion and the jet 
itself. This points to the fact that O stars with spherical 
fast winds would not be able in general to power some 
of the compact jets observed in the galaxy, implying that 
some special wind density profile should be given, likely 
produced by the presence of the compact object plus other 
factors like companion star rotation, etc (Paredes et al. 
2005J. For the launching radius, we have taken 4 i?sch to 
reach mildly jet Lorentz factors. Although our approach 
to estimate the energy balance between the jet and accre- 
tion is rough, this value for r\ is between the last stable 
orbit and the limits of the corona-like region, consistent 
with the present state of knowledge on this issue. 

6.1. Radio emission 

Previous models for leptonic emission from a magnetized 
compact jet predicted core dominated radio emission as 
observed in several galactic and extragalactic sources (for 
XRBs jet models, see, e.g. Markoff et al. 120011 for extra- 
galactic jet models, see, e.g. Ghisellini et al. .1989- O ur 
results are similar to those presented in these previous 
works, where radio emission that comes from the inner jet 
regions is strongly self-absorbed, becoming optically thin 
further down the jet. On the other hand, it is difficult 
to correctly model jet radio emission in some high-mass 
XRBs, like that presented by LS 5039 (Marti et al. IHMl 
Ribo et al. 2.999 ). A deeper discussion on this subject is 
presented in Paredes et al. H2005[) . We advance however 
that a z-dependence for the parameter rj, or also a less 
stringent restriction of £, could lead to a higher produc- 
tion of radio emission in optically thin regions of the jet. 

A correlation is seen between the luminosity in the ra- 
dio band and that at X-rays (Corbel et al. 120031 Gallo 
et al. 12003(1 that appears to be present in different sources 
with very different accretion rates and compact object 
mass values (Falcke et al.[2004). In the case of our model, 
when the source is jet dominated at X-ray wavelengths 
(cases A, B and D), the correlation is fulfilled. Otherwise, 
for corona-dominated sources (case C) the correlation can- 
not be reproduced, since we are not modeling the corona. 



Instead, we adopt a typical spectrum and a certain lumi- 
nosity in each case. 

7. Summary 

We have developed a detailed leptonic model for an MQ 
jet that can reproduce the emission observed from radio 
to gamma-ray energies, and makes precise predictions for 
very high energies. The basic assumptions of the model 
are a cold-matter dominated jet, with a magnetic field 
close to but below equipartition that is entangled with 
and dynamically dominated by jet cold matter. With the 
knowledge of the system parameters, given a simple stel- 
lar mass-loss density profile, and varying the jet to ad- 
vected matter ratio and the acceleration efficiency, a set 
of broad-band SEDs has been computed. Also, the opacity 
due to photon-photon interactions was taken into account 
to calculate the predicted SEDs. The absorption can sig- 
nificantly distort the production spectrum beyond 10 GeV 
mainly due to the effect of stellar photons in the case of 
massive companion stars. The opacity can be important 
even at lower energies when disk and corona radiation den- 
sities are high enough. The model shows that pair creation 
inside the jet could affect jet radiation. This aspect will be 
investigated accurately in future work. For systems where 
orbital eccentricity or other stellar mass-loss asymmetries 
are present, consistent predictions of the variability emis- 
sion pattern of the source can be obtained. 

New generation gamma-ray instruments, both 
satellite-borne or ground-based, like GLAST, AGILE, 
MAGIC or HESS can be used to test and constraint the 
high-energy predictions and assumptions of the model. 
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